High pressure phases of crystalhne tellurium. 
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A study of high pressure sohd Te was carried out at room temperature using Raman spectroscopy 
and Density Functional Theory (DFT) calculations. The analysis of the P-dependence of the experi- 
mental phonon spectrum reveals the occurrence of phase transitions at 4 GPa and 8 GPa confirming 
the high-pressure scenario recently proposed. The effects of the incommensurate lattice modulation 
on the vibrational properties of Te is discussed. DFT calculations agree with present and previous 
experimental data and show the metallization process at 4 GPa being due to the development of 
charge-bridges between atoms belonging to adjacent chains. A first-principles study of the stability 
of the 4 GPa phase is reported and discussed also in the light of the insurgence of lattice modulation. 

PACS numbers: 62.50.-p, 71.30.+h, 63.20.-e, 31.15.E-,61.44.Fw 



The remarkable advances in high-pressure (HP) meth- 
ods and techniques have surprisingly brought to hght the 
occurrence of incommensurately modulated (IM) crystal 
structures also in elemental systems under pressure 1- 
0]. The modulation, although preserving the long-range 
order, removes the lattice translational symmetry and in- 
troduces a new structural degree of freedom. Moreover 
in a variety of elemental systems (e.g. I, Br, P and VI 
group elements S, Se, Te) the picture is even more in- 
triguing since applied pressure induces a metallization 
transition preceding the onset of the IM phase, which 
suggests a connection with an electronic instability . 

Among the chalcogens the long-standing problem of 
the sequence and the nature of the HP phases of Te has 
been the subject of intense theoretical [1, and experi- 
mental [1, 0, investigations. Despite the great effort, 
a comprehensive picture is still far from being reached 
at least over the intermediate pressure range. It is in- 
deed well known that, under ambient conditions and on 
increasing the pressure up to 4 GPa, Te is a semiconduc- 
tor and its crystal structure is trigonal with spiral atom 
chains along the c axis (Te-I). The three atoms in the 
unit cell are linked via covalent-like bonds to the nearest 
neighbors along the chain and, via Van der Waals in- 
teractions, to the second neighbors lying along adjacent 



chains [ll|. The very HP regime (i.e. above 27 GPa or 29 
GPa according with Ref. [l2| or Ref. @ respectively) is es- 
tablished, as well: Te is a metal with a highly symmetric 
bcc structure (Te-IV). 

Apart from the metallization transition that, since 
long, it is known to occur at 4 GPa [13] together with 
a structural change, there is an intense debate on the 
HP metallic structural phases between the Te-I and the 
Te-IV. i.e. over the 4-27(29) GPa pressure range. In 
particular the succession of phases over this intermedi- 



ate pressure range accepted up to few years ago was: 
a puckered layer monoclinic structure for P > 4 GPa 
which transforms into an orthorhombic one at 6.8 GPa 
and finally a /3-Po-type above 11 GPa and below the Te- 



IV pressure threshold |lll - ll3j . Recent X-ray diffraction 
studies d, [l3| have drawn a rather different HP phase 
transition scenario, claiming a structural transition to 
triclinic (Te-II) at 4.0 GPa, followed by the onset of an 
IM monoclinic phase (Te-III) at 4.5 GPa which is finally 
transformed into the bcc Te-IV phase above 29 GPa. It is 
worth to notice that since Te-I is observed up to 4.5 GPa 
and Te-II up to 8 GPa, we have large regions of coexist- 
ing phases and the above reported phase boundaries are 
intended as the pressure threshold at which the higher 
pressure phase first appears on compressing the lattice. 

In the present paper a combination of experiments and 
ab initio calculations is used to study structural, dynamic 
and electronic properties of HP solid Te in order to solve 
the controversy about the sequence of crystallographic 
structures, the real occurrence of IM phases, and to shed 
light on the microscopic metallization mechanism. A HP 
Raman study of Te (0 — 15 GPa) is reported together 
with the results of density functional theory (DFT) cal- 
culations which well reproduce the experiment and allow 
to gain a further insight of the microscopic modifications 
induced by pressure. 

A high purity Te sample (99.99% by Aldrich) was used 
without further treatment. A diamond anvil cell coupled 
to a microRaman spectrometer (Horiba Jobin Yvon) was 
employed. The ruby fluorescence technique was applied 
for pressure calibration and two pressure-transmitting 
media (NaCl and methanol/ethanol 1:4 mixture) were 
used to verify their possible effects on the high pressure 
measurements. Moreover, exploiting the high spatial res- 
olution of our experimental setup, Raman spectra were 
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collected at each pressure from different points over the 
surface of the samples (about 50x50 /xm^) to test the 
presence of pressure gradients. Further experimental de- 
tails can be found in Ref. [isj, here we just notice that 
great care was taken in the optical alignment [31 to col- 
lect reliable spectra above 80 cm~^ and that we did not 
find any evidence of either reactions between the me- 
dia and the sample or detectable non-hydrostaticity ef- 
fects. No laser-induced sample heating has been observed 
as well[l3l. Normalized and background subtracted (see 
below) Raman spectra of Te collected using NaCl and 
methanol/ethanol are shown in Fig. [l]at selected P val- 
ues. According to the standard group theory (GT), the 
irreducible representation of the optical modes of Te-I is 
Ai -I- A2 -f- 2 E. The Raman active modes are: the Ai 
mode (breathing in the ab plane) and the two doubly 
degenerate E modes (E(l) a- and 6-axis rotation , E(2) 
asymmetric stretching mainly along c-axis) . These latter 
are also Infrared active since the crystal lacks a center of 
inversion. Measured spectra were fitted by using a stan- 
dard model curve [l8| given by the sum of three phonon 
contributions (each described by a damped harmonic os- 
cillator), and an electronic background. A linear baseline 
was also included in the fitting curve. In Fig. [T]d the best- 
fit curve and the phonon assignment are also shown for 
the lowest pressure spectrum. The spectra shown in Fig. 
[T^ and Fig. [T}d are background subtracted and normal- 
ized to the phonon integrated intensities. 

The Raman spectrum shows a remarkable pressure de- 
pendence albeit the three-peaks structure is basically pre- 
served over the whole pressure range. We notice that the 
persistence of a Raman signal well above 11 GPa dis- 
agrees with the old phase transition scheme 11- 13| since 
no Raman-active are expected for the /3-Po-type phase. 
On the other hand, an unmodulated monoclinic phase 
(i.e. an unmodulated Te-III structure) has no Raman- 
active mode as well. The observation of a Raman sig- 
nal is thus only coherent with the new transition scheme 
where the Te-III structure appears to be IM Q. In this 
case, optical scattering processes are allowed not only 
from phonons at F but, in the sinusoidal approximation, 
also at fc = ±<7i (with qt the modulation wave vector). 
Using the superspace GT (equations and formalism are 
from Ref. [IB]): the calculated character of the represen- 



tation is xi^j 



, /) = (6, 2, 0, 0) where E, 2y, uiy de- 



note the point group operations of the three-dimensional 
space-group elements. This representation is reducible 
into F = 2 Ag -I- 1 Bg 4- 2 A„ -1- 1 B„, where three Ra- 
man active modes (two Ag modes and the Bg mode) are 
expected. This first experimental evidence thus provides 
not only a clear support to the recently proposed transi- 
tion scheme but can be seen as a marker of the onset of 
an IM phase. 

Looking at Fig. [TJ;, three different pressure regions can 
be identified by the abrupt variations of the pressure de- 
pendence of the phonon peak frequencies. The borders 
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FIG. 1. (color online) Te Raman spectra at RT and selected 
P, in ethanol- methanol (a), and NaCl (b). Best fit curve and 
phonon contributions for the P=1.8 GPa spectrum are also 
shown in (b). P-dependence of phonon frequencies (c) and 
linewidths (e). Vertical lines in (c) and (e) mark the phase 
transitions pressures according to ref. \XM- The intensity 
Iratio = /[£;(2)]/(7[_B(l)] -l-/[^i]) (see text) vs. pressure in 
(d). Open (close) symbols in (c), (d) and (e) refer to measure- 
ments using ethanol-methanol mixture (NaCl). Color code 
and atomic displacements for the Raman active modes in Te- 
1(f). 



among the regions actually correspond to the pressure 
thresholds of the new transition scheme [l3| • Over the 
0-4 GPa pressure range, the phonon frequency of the Ai 
mode shows a linear softening (~ 10 cm~^), whereas the 
frequencies of the two E modes weakly increase. Since the 
Ai mode is a chain breathing and, over the same pres- 
sure range, a contraction of the b lattice parameter and 



an almost constant c 11| were observed, the above find- 



ings suggest a progressive weakening of intrachain bonds 
in favor of interchain atomic interactions. The slope dis- 
continuity of the pressure dependence of the E(2) mode 
whose frequency starts to decrease at 4 GPa, marks the 
emergence of a new crystal phase albeit, within the 4-7 
GPa, the frequency of the Ai phonon is still decreasing 
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and that of the the E(l) mode keeps almost constant. 
The narrow pressure range in which Te-I and Te-II coex- 
ist (4-4.5 GPa) makes it impossible to determine whether 
the discontinuity can be ascribed to the Te-I/Te-II tran- 
sition or to the onset of the Te-III phase. According to 
Ref. [1, that reports on the Te-II/Te-III coexistence 
within the 4.5-8 GPa pressure range, the observed soft- 
ening of the E(2) mode can be ascribed to the structural 
instability of the Te-II phase which is progressively con- 
verted into the rather similar Te-III phase We 
notice that larger values of the phonon linewidths shown 
in Fig. [1^ can be found within the coexistence region 
(4.5-8 GPa) as a consequence of the lattice disorder. On 
further increasing the pressure above ~ 8 GPa all the 
phonon modes show a rather regular frequency harden- 
ing, indicating that a single, stable phase (Te-III) has 
finally established. 

The presence of modulation implies a one-to-one corre- 
spondence between eigenvalues of the dynamical matrix 
at generic k and those at fc ± More precisely the 

phonons at P involved in the extra scattering processes 
arising from the modulation, have the same energies as 
the phonons at the point qi involved in the usual inelastic 
scattering processes. [3| This allows an unusual compari- 
son between the results of a Raman experiment and those 
of an inelastic x-ray (neutron) scattering experiment. In- 
deed, the frequency of the Raman mode at higher en- 
ergy, observed at 8.9 GPa (Te-III phase, see Fig. [Tt), 
is the same within the experimental uncertainties as the 
frequency of the longitudinal mode at k = qi observed 
by the inelastic x-ray scattering experiment reported in 
Ref. ,7]] . This finding thus provide further evidence of the 
the onset of an IM phase at HP. 



GPa 




FIG. 2. (color online) Valence electron density at and 5 
GPa in the Te-I phase over the plane containing 3 atoms along 
the chain and the closest atom in the adjacent chain. 

The weakening of intrachain bonds can be ascribed 
to the onset of a pressure induced charge transfer pro- 
cess, as also suggested by the remarkable reduction of 
the intensity (peak integrated area) of the Ai mode, 
I[Ai] [13 . This effect is shown in Fig. [TJi where 
the pressure dependence of the intensity ratio Iratio = 



I[E{2)]/{I[Ai]+ I[E{1)]) is shown. Neglecting the smah 
contribution to Iratio due to the small intensity (area) 
of the E(l) mode, the growth of Iratio indicates a strong 
decrease in the intensity of the Ai peak before the metal- 
lization pressure. Moreover two different regimes, below 
and above the metallization pressure (4 GPa), can be 
clearly identified. 

To gain further insight into the microscopic charge- 
transfer process, pseudopotential DFT calculation us- 
ing the Quantum-Espresso code [2l| have been carried 
out. A full structural relaxation under the influence of 
Hellmann-Feynman forces and stresses was carried out 
for Te-I at 0, 2.5 and 5 GPa using the variable cell shape 
method (221]. At each pressure the electronic density of 
states and bands as well as the frequencies and the inten- 
sities of normal vibrational modes were obtained by the 
perturbation approach of DFT [2^ . We found the closure 
of the energy gap at 5 GPa and an overall good agree- 
ment between the calculated data and the HP experi- 
mental structural data ll|. The intrachain bond length 



was slightly overestimated because of the typical error of 
the generalized gradient approximations in the determi- 
nation of the cohesive force in strongly anisotropic sys- 
tems. This reflects in an underestimate of the calculated 
phonon frequencies (by less than 10%) which, however, 
together with the intensity ratio Iratio defined above, 
shows the same pressure dependence as the experimental 
data. Moreover, the calculations support the evidence 
for a progressive reduction of the angle between the E(2) 
mode eigenvectors and the c~axis on applying pressure 
(from 12° at P=0 to 2.5° at P=5). The good agreement 
of the results of DFT calculations with the present and 
the previous experimental data makes us confident also 
about the other possible outputs of the code. In par- 
ticular our calculations provided also the maps of the 
valence electron (5S^5P''^) densities shown in Fig. [5] at 
ambient pressure and at 5 GPa . The onset of a bridge 
of charge density linking the atoms inside the chain and 
the first neighbor atom in the adjacent chain is pretty 
clear at P=5 GPa. The pressure induced metallization 
transition is therefore the result of the simultaneous lat- 
tice symmetrization and charge-transfer processes from 
intrachain to interchain regions which allows the mobil- 
ity of the valence electron. 

We also exploited the DFT approach to investigate the 
structural Te-I/Te-II transition although it has proved to 
be a rather complex task. In Te-I most of the crystallo- 
graphic constants are fixed by the hexagonal geometry, 
nevertheless the DFT underestimates the anisotropy of 
the system leading to a slight overestimation of the ratio 
of chain radius to interchain distance. In the Te-II case, 
where all the lattice parameters and the atom positions 
are left free to vary during the cell relaxation, this trend 
is strongly emphasized, so that, with relaxing the struc- 
ture, the experimental zig-zag puckering is entirely lost. 
For this reason, the dynamic properties for the Te-II re- 
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laxed cells in the pressure range between 3.5 and 7 GPa 
differ qualitatively from the experimental ones, even if 
the calculated enthalpy of the Te-I phase becomes higher 
than that of Te-II phase above 5.5 GPa, in agreement 
with the occurrence of the phase transition (see Fig. [3]) . 
Thus, in the present case, the dynamical properties were 
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FIG. 3. The difference AH = Hre-ii - Htb-i between the 
enthalpy values of the Te-II and Te-I cells relaxed at different 
pressures using the variable cell shape method. Dotted lines 
mark the crossover. 



to an excessive sensitivity of the DFT theory to predict 
the incommensurate phase, or in other words it seems to 
be related to an overestimate of the width of the sad- 
dle point of the energy landscape in the direction of the 
modulation. This is confirmed by the fact that there is a 
neighborhood of points around q^c for which the trans- 
verse acoustic phonon has negative energies [6j . 

In summary the Raman experiment supports the re- 
cently proposed HP phase transition scenario Q and the 
appearance of an incommensurate modulation is clearly 
demonstrated at intermediate pressure. Our DFT cal- 
culations well reproduce the present and the previous 
HP experimental data allowing for an accurate study 
of the microscopic pressure induced modifications. In 
agreement with the indications of the Raman experiment, 
DFT calculations show that the metallization transition 
is driven by an intra- to inter-chain charge transfer which 
develops charge density bridge among adjacent chains. 
Finally the combined use of spectroscopic and ah ini- 
tio methods has proved to be a useful approach to the 
complexity of systems with incommensurately modulated 
structure. 



calculated optimizing the atomic positions while keeping 
fixed the lattice parameters to the experimental values 
obtained for the Te-II phase at 4.5 GPa (iJ]. For sake 
of comparison the Raman frequencies have been recal- 
culated by optimizing only the atomic positions keeping 
constant the lattice parameters to the experimental val- 
ues observed for the Te-I phase at 4.0 GPa. The results 
for the two phases are reported in Tab. H] 



TABLE I. Calculated Raman phonon frequencies for Te-I and 
Te-II optimizing the atom positions while keeping constant 
the lattice parameters to the experimental values. 



Te-I at 4 GPa 


Te-II at 4.5 GPa 


Mode 


u (cm ^) 


Mode 


u (cm ^) 


E(l) 


78 


As(l) 


38 


Ai 


92 


A«(2) 


94 


E(2) 


122 


As (3) 


125 



We notice that only the low-frequency phonon shows a 
remarkable discontinuity at the transition whereas the 
differences between the values calculated in the two 
phases for the high-frequency modes are within the typ- 
ical DFT calculations accuracy. In agreement with the 
experimental data our calculations thus predict a con- 
tinuous evolution of the frequencies across the Te-I Te- 
II transition for the high-energy modes. As to the es- 
timate of the frequency of the low-energy mode further 
discussion is required. This vibrational mode is the mode 
mostly affected by the IM since the the direction of the 
incommensurate modulation is rather close to the direc- 
tion along with the atom displacements occur The un- 
derestimate of the mode frequency seems to be related 
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